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ABSTRACT 

We investigate the formation of disk-bulge-halo systems by including bulges in 
the Fall & Efstathiou theory of disk formation. This allows an investigation of bulge 
dominated disk galaxies, such as SOs and disky ellipticals. These latter systems, which 
consist of an elliptical spheroid with an embedded disk with a scale-length of typically 
a few hundred parsecs, seem to form a smooth sequence with spirals and SOs towards 
lower disk-to-bulge ratio. The aim of this paper is to examine whether spirals, SOs, and 
disky elllipticals all can be incorporated in one simple galaxy formation scenario. We 
investigate an inside-out formation scenario in which subsequent layers of gas cool and 
form stars inside a virialized dark halo. The inner, low angular momentum material is 
assumed to form the bulge. Stability arguments are used to suggest that this bulge 
formation is a self-regulating process in which the bulge grows until it is massive 
enough to allow the remaining gas to form a stable disk component. We assume that 
the baryons that build the disk do not loose their specific angular momentum, and 
we search for the parameters and physical processes that determine the disk-to-bulge 
ratio, and therewith explain to a large extent the origin of the Hubble sequence. The 
spread in halo angular momenta coupled with a spread in the formation redshifts can 
explain the observed spread in disk properties and disk-to-bulge ratios from spirals to 
SOs. If galaxy formation is efficient, and all available baryons are transformed into the 
disk-bulge system, cosmologies with ^ 0.3 can be excluded, since stable spiral disks 
would not be allowed to form. However, if we assume that the efficiency with which 
galaxies form depends on the formation redshift, as suggested by the small amount of 
scatter in the observed Tully-Fisher relation, and we assume that the probability for a 
certain baryon to ultimately end up in the disk or bulge is independent of its specific 
angular momentum, spirals are allowed to form, but only at small formation redshifts 
[z <^ 1). At higher formation redshifts, stability arguments suggest the formation of 
systems with smaller disk-to-bulge ratios, such as SOs. Since density perturbations in 
clusters will generally collapse earlier than those in the field, this scenario naturally 
predicts a density-morphology relation, the amplitude of which depends on the baryon 
fraction of the Universe. Disky ellipticals are too compact to be incorporated in this 
scenario, and they thus do not form a continuous sequence with spirals and SOs, 
at least not in the sense of the galaxy formation scenario envisioned in this paper. 
Alternative formation scenarios for the disky ellipticals, such as gas-rich mergers or an 
internal mass loss origin for the embedded disks, are much more viable. 

Subject headings: galaxies: formation — galaxies: general — galaxies: halos - 
instabilities 
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1. Introduction 

One of the most compelling puzzles in present day astronomy is the question how galaxies 
formed. In particular, we need to understand the wide variety of sizes, masses and morphologies 
of galaxies observed, as well as their dynamics, ages and metallicities. In addition, we need 
to understand the origin of the scaling relations, such as the fundamental plane relations for 
ellipticals and the Tully-Fisher relation for spirals, as well as the so-called density- morphology 
relation (Dressier 1980). The latter shows that the more compact galaxies, such as ellipticals and 
SOs, are preferentially found in overdense regions such as galaxy clusters. 

The classical way of depicting the different morphologies of galaxies is by means of the 
Hubble diagram, which reveals the gross distinction of three sorts of galaxies: spirals, SOs, and 
ellipticals. Several studies over the past years have shown that elliptical galaxies can be roughly 
divided in two subclasses: the rotation supported, low luminosity disky ellipticals and the pressure 
supported, bright, boxy ellipticals (e.g., Bender 1988; Bender et al. 1989; Capaccioli, Caon & 
Rampazzo 1990). This dichotomy has recently been strengthened by properties observed in 
their central regions (e.g., Nieto, Bender & Surma 1991; Jaffe et al. 1994; Ferrarese et al. 1994; 
Gebhardt et al. 1996; Faber et al. 1997). Based on this dichotomy amongst elliptical galaxies, it 
has been suggested that the classical Hubble diagram should be revised (Kormendy & Bender 
1996), to use a more physical subdivision of the class of the ellipticals, rather than the apparent 
flattening. Kormendy & Bender proposed to use the velocity anisotropy which is well measured 
by the isophote shapes (Bender 1988; Bender et al. 1989). 

The main morphological parameter that sets the classification of galaxies in the (revised) 
Hubble diagram is the disk-to-bulge ratio D/B. Understanding the origin of the Hubble sequence 
is thus intimately related to understanding the parameters and processes that determine the ratio 
between the masses of disk and bulge. The diskiness of the low luminosity ellipticals is generally 
interpreted as due to an embedded disk. These disks, which we term 'embedded disks' in the 
following, are smaller and brighter than disks in SOs and spirals (Scorza & Bender 1995), and it 
has been emphasized many times that the disky ellipticals build a continuous sequence with SO 
galaxies (Capaccioli 1987; Carter 1987; Bender 1988, 1990; van den Bergh 1990: Capaccioli et 
al. 1990; Rix & White 1990; Capaccioli & Caon 1992; Bender et al. 1993; Scorza & Bender 1995, 
1996; Scorza & van den Bosch 1998). This continuity suggests similar formation histories, whereby 
one or several parameters of the proto-galaxy vary smoothly. 

Disk dominated systems such as spirals are believed to have formed by cooling of the baryonic 
matter inside a virialized dark halo. As the gas cools, its specific angular momentum is conserved, 
and the amount of angular momentum of the dark halo thus determines the size of the disk (Fall 
& Efstathiou 1980; see Section 2 for more details on this disk- formation scenario). The formation 
of bulge dominated disk systems is far less clear. The aim of this paper is to investigate to what 
extent the disky ellipticals and SOs can be incorporated in the Fall & Efstathiou theory for the 
formation of galactic disks, by incorporating a simple picture for the formation of the bulge. 

We envision an inside-out formation scenario for the bulge. It is assumed that the bulge 
forms out of the low-angular momentum material in the halo, which cools and tries to settle into 
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a small, compact disk. Such disks are however unstable, and we assume here that this instability, 
coupled with the continuous supply of new layers of baryonic matter that cool and collapse, forms 
the bulge. This inside-out bulge formation is self-regulated in that the bulge grows until it is 
massive enough to allow the remaining gas to form a stable disk component. We do not describe 
the bulge formation in any detail but merely use empirical relations of the characteristic structural 
parameters of bulges and ellipticals to describe the end result of the formation process as a 
realistic galaxy. We use this simple formation scenario to investigate the predicted disk-to-bulge 
ratios and disk scale-lengths as a function of the halo angular momentum, and as a function of 
formation redshift and cosmology. The main focus of this paper is to investigate whether this 
inside-out formation scenario can account for two orders of magnitude variation in disk-to-bulge 
ratio, required in order to incorporate a wide variety of disk-bulge systems: from late-type spirals 
with D/B £ 10 to SOs (D/B ~ 0.1) to disky ellipticals (D/B ~ 0.1). 

This paper is organized as follows. In Section 2 we describe the formation scenario, the galaxy 
formation efficiency, and the stability criterion for our disk-bulge-halo models. In Section 3 we 
use these models to investigate the position of different sorts of disks in the parameter space of 
disk central surface brightness versus disk scale-length. In Section 4 we discuss to what extent 
SOs and disky ellipticals can be incorporated in this formation scenario, and what parameters 
may be responsible for the origin of (the major part of) the Hubble sequence. In Section 5 we 
briefly discuss alternative formation scenarios for disky ellipticals. Our results are summarized 
and discussed in Section 6. 



2. The formation scenario 

A remarkably successful model for the formation of galaxies is the White & Rees (1978) 
theory, wherein galaxies form through the hierarchical clustering of dark matter and subsequent 
dissipational settling of gaseous matter within the dark halo cores. Coupled with the notion of 
angular momentum gain by tidal torques induced by nearby proto-galaxies (Hoyle 1953; Peebles 
1969), the White & Rees theory provides the background for a model for the formation of galactic 
disks (Fall & Efstathiou 1980; Faber 1982; Gunn 1982; Fall 1983; van der Kruit 1987). In this 
model, disks form through the collapse of a uniformly rotating proto-galaxy. Owing to the 
non-dissipative character of the dark matter, its collapse halts when the system virializes. At 
regions where the baryonic density is high enough, the baryons can cool and decouple from the 
dark matter. The angular momentum of the baryons, acquired from tidal torques from nearby 
proto-galaxies, prevents the collapse from proceeding all the way to the center, and causes the 
baryons to settle in a rapidly rotating disk. This simple picture of disk formation yields realistic 
disk sizes and surface brightness profiles provided that there is little angular momentum transfer 
from the baryonic matter to the dark halo, and that the dark halo has a density profile such that 
the resulting galaxy has a flat rotation curve. Fall (1983) extended the original Fall & Efstathiou 
disk formation scenario by placing it in a cosmological context, and by deriving Freeman's law 
(Freeman 1970) and the Tully-Fisher relation. Recently, two papers reexamined this disk formation 
model and included the effects of halo-contraction due to the accumulation of baryonic matter in 
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the center of the galaxy: Dalcanton, Spergel & Summers (1997) showed that the disk formation 
model outlined above can not only explain the properties of normal, high surface brightness 
(HSB) disks, but also those of the class of low surface brightness (LSB) disks (see also Jimenez et 
al. 1998); Mo, Mao & White (1998) compared the outcome of disks in different cosmologies, and 
used this to investigate the effects of disk evolution and the predicted population of Lya absorbers 
in QSO spectra. The assumptions we make are similar to those made by Dalcanton et al. (1997) 
and Mo et al. (1998), but we extend upon their work by including bulges. We investigate how the 
disky ellipticals and SOs may be fitted into this particular formation scenario, by comparing the 
differences in properties of the proto-galaxies that result in spirals, SOs, and disky ellipticals. 



2.1. Inside-out bulge formation 

The turn-around, virialization, and subsequent dissipational settling of the baryonic matter of 
a proto-galaxy is an inside-out process. First the innermost shells virialize and heat its baryonic 
material to the virial temperature. Because of the relatively high density and low angular 
momentum of this gas, it will rapidly cool and (try to) settle in a very small, compact disk. 
However, self-gravitating disks are violently unstable and form a bar (e.g., Hohl 1971; Ostriker & 
Peebles 1973). Meanwhile the next shell of the proto-galaxy cools and tries to settle into a disk 
structure at a radius determined by its angular momentum. The resulting structure in the center 
of the dark halo consequently becomes more and more unstable. 

Bars are efficient in transporting gas inwards (e.g., Friedli & Martinet 1993; Wada & Habe 
1995), and can cause vertical heating by means of a collective bending instability (e.g., Combes 
et al. 1990; Pfenniger 1984; Pfenniger & Friedli 1991; Raha et al. 1991). Both these process lead 
ultimately to the dissolution of the bar; first the bar takes a hotter and triaxial shape, but is later 
on transformed in a spheroidal bulge component (e.g., Combes et al. 1990; Pfenniger &: Norman 
1990; Pfenniger & Friedli 1991; Pfenniger 1993; Friedli 1994; Martinet 1995; Norman, Sellwood & 
Hasan 1996). Here we assume that this bar instability and subsequent dissolution, coupled with 
the continuous supply of new layers of baryonic matter that cool and collapse, forms the bulge. 
The numerical simulations performed by the numerous studies listed above, have been able to 
show that bars can build small and close to exponential bulges (see also Courteau, de Jong & 
Broeils 1996 for observational evidence for bulge formation out of disks). However, no simulations 
have been able to use bar instabilities to produce very large bulges or even ellipticals. On the other 
hand, all these simulations started from a normal, marginally unstable disk, whereas the process 
we are describing starts from a highly unstable, compact disk, includes the continuous supply of 
baryonic matter from prolonged cooling flows, and is likely to have additional heating sources due 
to the star formation and feedback processes which are expected to take place in this phase of the 
bulge formation. Detailed numerical simulations that take all these processes into account are 
required in order to investigate whether this can indeed lead to massive bulges as assumed here. 
Finally, we emphasize that the high densities of the inner shells yield cooling time-scales that may 
be significantly shorter than the dynamical time-scale. This can add to the inside-out formation 
of a bulge component, as suggested by Kepner (1997). 
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We assume that the bulge formation process discussed above continuous until the bulge 
has become massive enough such that the subsequent layers of baryonic material can cool and 
form a disk which is stable against bar-formation. The process of disk-bulge formation is thus a 
self-regulating one in that the bulge grows until it is massive enough to sustain the remaining 
gas in the form of a stable disk. Since the stability is directly related to the amount of angular 
momentum of the gas (see Section 2.4), we expect a clear correlation between the disk-to-bulge 
ratio and the angular momentum of the dark halo out of which the galaxy is assembled. We will 
show that this correlation is present, and that it nicely follows the marginal stability curve as 
expected from the self-regulating mechanism proposed here. 



2.2. A disk-bulge-halo model 

The galaxies under investigation in this paper, varying from spirals to disky ellipticals, are 
assumed to consist of a disk plus bulge embedded in a dark halo. Virialized dark halos can be 
quantified by a mass M, a radius r2oo, and an angular momentum J (originating from cosmological 
torques). Spherical collapse models (e.g., Gunn & Gott 1972) show that virialized halos have an 
over-density of approximately 200. In a recent series of papers Navarro, Prenk & White (1995; 
1996; 1997, hereafter NFW) showed that virialized dark halos have universal equilibrium density 
profiles, independent of mass or cosmology (but see Moore et al. 1997), which can be well fit by 

p(r) = Pent . ■ r -T2 > (1) 

(r/r s )(l +r/r s y 

where 

_ 200 c 3 
3 ln(l + c) - c/(l + c) ' 1 ' 

with c = r2oo/r s the concentration parameter, r s a scale radius, and p cr \t the critical density for 
closure. The radius r2oo is defined as the limiting radius of a virialized halo, and corresponds to 
the radius within which the mean density is 200p cr i t . Throughout we only consider spherical halos. 
Rather than specifying the halo by its mass M, it is customary to specify it by its circular velocity 
V200 = v / G , M/r 2 oo- Given the mass M, the formation redshift £f orm , and the specific cosmology, 
the concentration parameter c can be calculated using the procedure outlined in Appendix A of 
Navarro, Frenk & White (1997). In this paper we compare two different cosmologies: a standard 
cold dark matter cosmology (SCDM) with a total matter density of 0,q = 1.0, a Hubble constant 
of Hq = lOO/i kms _1 Mpc _1 with h = 0.5, and with a rms linear over-density at z = in spheres 
of radius 8/i _1 Mpc of dg = 0.6. This latter parameter normalizes the CDM power spectrum. 
The second cosmology discussed here is an open cold dark matter model (OCDM) with Qq = 0.3, 
h = 0.6, and ag = 0.68. Both models have a zero cosmological constant and the parameters 
are typical of those favored by large-scale structure constraints. The baryon density is set to be 
^bar = 0.0125 /i -2 , in agreement with the nucleosynthesis constraints (e.g., Walker et al. 1991). 
The baryonic mass fraction is then given by f^bar/^o- 



We assume the final system to have a disk-to-bulge ratio t = M^/Mb, such that the masses of 
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disk, M d , and bulge, Mb, can be written as 

Md = e af -^r t iff M = m ^ ( 3 ) 

and 

Mb = e 9f YT~ t iff M = mbM - (4) 

Here e g f is the galaxy formation efficiency, expressed as the fraction of the available baryons that 
ultimately end up in the galaxy (see Section 2.3). As mentioned in Section 2.1, we assume that 
the bulge is formed out of the low-angular momentum material. Furthermore, we assume that 
the material that builds the disk does not loose its specific angular momentum. Under these 
conditions we have that J d = j d J, where 

H = f(t,e g f)[m d + m b ]. (5) 

The function f(t, e g f) depends on the disk-to-bulge ratio and on the detailed physics that describe 
how the fraction e g f of the baryonic material ends up in the disk-bulge system (see Section 2.3). 
This function is derived in the Appendix. In the following we parameterize J by the dimensionless 
spin parameter A, defined by 

X = GM^- (6) 
Here E is the halo's total energy. Several studies, both analytical and numerical, have shown that 
the distribution of spin angular momentum of collapsed dark matter halos is well approximated 
by a log normal distribution: 

1 ( ln 2 (A/A)\dA 

(e.g., Barnes & Efstathiou 1987; Ryden 1988; Cole & Lacey 1996; Warren et al. 1992). We consider 
a distribution that peaks around A = 0.05 with a\ = 0.7. These values are in good agreement with 
the iV-body results of Warren et al. (1992). 

We assume the disk to be an infinitesimally thin exponential disk, with surface brightness 

£ d (12) = Z dfi exp{-R/R d ). (8) 

and total luminosity 

L d = 27r£ di0 i4 (9) 

In order to convert the disk parameters to an observers frame we need to assume a mass-to-light 
ratio T d . Bottema (1997) found T B = (1.79 ±0.48) for h = 0.75 for HSB disks. Using B — V » 0.8 
(e.g., de Jong 1996b) and (B — V) & = 0.65 this yields = 2.1h in the U-band. 

Although bulges and elliptical galaxies are reasonably well fit by a r 1//4 'de Vaucouleurs' law, 
it has recently been shown that a better fit is obtained by using a generalized r l l n law. This 
fitting function, introduced by Sersic (1968), has the form 

l/n" 



E b (r) = S bj0 exp 



-«(-) 



(10) 
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where £& o is the central surface brightness, r e is the effective radius, and n is the exponential 
variable. For n = 1 one has a pure exponential, while for n = 4 one recovers the de Vaucouleurs 
profile. The parameter o is determined by the requirement that r e is the radius encircling half of 
the total luminosity and is very well approximated by 

a = 2.0n - 0.324, (11) 

(Ciotti 1991). The total luminosity of a system with a r l l n luminosity profile (assuming sphericity) 
is given by 

nT(2n) „ o , 
L b = -^^^bfirl 12 

where V is the Gamma function. It was found that the best fitting exponent n is correlated with 
the total luminosity of the spheroid: Andredakis, Peletier &: Balcells (1995) used the r 1//n law to fit 
the luminosity profiles of bulges of spiral galaxies, and found the effective radius r e to be related 
to the total luminosity of the bulge by the empirical relation 

M B = -19.75 -2.8 log(r e ), (13) 

In addition, Caon, Capaccioli & D'Onofrio (1993) found that the exponent n is correlated with 
the effective radius through the empirical relation 

log(n) = 0.28 + 0.52 log(r e ). (14) 

Although this empirical relation has only been shown to fit early-type galaxies, we assume here 
that it is also valid for bulges. We take n to be limited between n = 1 (exponential) and n = 4 (de 
Vaucouleurs). None of our results are significantly influenced by this assumption however. The 
main parameter for the bulge is its total mass; changes in its actual density distribution are only a 
second order effect. The relations (|l0|) - (|i~4]) allow us to describe the bulges by a single parameter; 
once the total luminosity of the bulge is known, we can use the empirical relations to determine 
all other bulge parameters! Van der Marel (1991) found that ellipticals have T# = 6.64/i (Johnson 
R). Using V - Rj « 0.68 (e.g., van der Marel 1991) and (V - R j) Q = 0.52 this yields T b = 7.7h in 
the U-band. For our combined disk-bulge model we assume that the mass-to-light ratio is constant 
over the entire galaxy with a value equal to the luminosity weighted average of the mass-to-light 
ratios of the disk and bulge. Alternatively, we could have assigned the bulge and disk separate 
mass-to-light ratios of and T;, respectively, but the mass-to-light ratio of the disks in SOs and 
disky ellipticals is clearly larger than for spiral discs. Assuming a constant, luminosity weighted 
ratio is thus more realistic. 

Once the baryons inside the virialized dark halo start to cool, this induces a contraction of 
the inner region of the dark halo. We follow Blumenthal et al. (1986) and Flores et al. (1993) 
and assume that the halo responds adiabatically to the slow assembly of the baryonic matter in 
disk and bulge, and that it remains spherical (see also Dalcanton et al. 1997 and Mo et al. 1998). 
Under these assumptions the angular action of the dark halo material is an adiabatic invariant, 
such that a particle initially at a radius r, ends up at mean radius r according to 



r i M i (r i )=rM f {r). 



(15) 
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Here Mj-(r) is the initial mass distribution (the NFW profile) given by 

ln(l + cx) — cx/(l + cx) 



Mi(r) = M 



ln(l + c) - c/(l + c) 



with x = r/r2oo, and M/(r) is the final mass distribution given by 



Here 



and 



M f (r) = M d (r) + M b (r) + (1 - m d - m 6 )M i (r i ). 
M d (r) = M d 



M 6 (r) 



2a 2 » +1 
7rn 2 r(2n) 



1 - (l + — ) eM-r/Rd) 



r/r e oo 

M b J dyy 2 J x 13 ' 1 exp[-ax p ] 



dx 



V^ 2 



If 



with = 1/n. 

Mo et al. (1998) have shown that the scale- length of the disk can be written as 

^ = 4j(— )Ar 20 o/ c - 1/2 //, 



where 



fR 2 



^200 



a parameter that indicates the amount of self-gravity of the disk, and 

_ c 1 - 1/(1 + cf - 2 ln(l + c)/(l + c) 
/c " 2 [c/(l + c) - ln(l + c)] 2 

The circular velocity V c (r) is given by 

V c \r) = V c 2 d (r) + V 2 b (r) + y c 2 DA/ W' 

with 



GM b {r) 



V 2 d (r) as given by equation [2.196] in Binney & Tremaine (1987), and 



V, 



2 , ,_ G[Mf(r)-M d (r)-M b {r)} 

c,DMV) — ~ 



(16) 

(17) 
(18) 

(19) 



(20) 
(21) 

(22) 

(23) 
(24) 

(25) 



We can solve this set of equations in an iterative way, following the procedure used by Mo 
et al. (1998). For a given M, A, Zf orm and t we first calculate the concentration parameter c of 
the dark halo before collapse of the baryons. We then calculate the parameters of the bulge using 
equations (|i~T| ) - fli~4]). Next we start with a guess for by setting /r = 1, and we calculate the 
disk parameters R d and Eo,d using equations ( |20| ) and @ respectively. Subsequently we solve for 



rj as a function of r using equations ( p^5[ ) - (|19j), and thus obtain Mf(r). We then calculate 



and substitute this in equation ( |21] ) to obtain a new value for /r. Convergence is rapidly achieved: 
3 to 5 iterations yield an accuracy of R d already better than one percent. 
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2.3. Estimating the galaxy formation efficiency 

Since we have denned the mass of the dark halo such that the average over-density of the halo 
is 200 times the critical density for closure at the virialization redshift z, we have that 

M = V ™, . (26) 
10GH(z) v ; 

This halo mass can be related to the disk luminosity (assuming no bulge, i.e., t = oo), by 

e 9 /(^bar/^o) ' ^ 

(cf. equation ||). Combining the above equations yields 

_ 6g/(^bar/^o) T/3 / 9 o>, 

Ld ~ WGT d H(z) V20 °- (28) 

This is similar to the Tully-Fisher relation: L = A(y c /220)" 1 , whereby the zero-point, A, is 
inversely proportional to H(z). However, current data suggest no zero-point evolution (Vogt et 
al. 1996, 1997; Mao, Mo & White 1997). In other words, if we assume a significant scatter in the 
formation redshifts of spiral galaxies (some massive spirals are already in place at z ~ 1, see Vogt 
et al. 1996), the dependence of equation (pq) would yield a scatter in the Tully-Fisher relation 
much larger than observed: if spirals form between z = and z = 1 (and they are not significantly 
modified since their formation) a scatter of ~ 1.1 mag is predicted for an Einstein-de Sitter 
Universe, much larger than the observed scatter of <^ 0.3 mag. This suggests that the redshift 
dependence of the Tully-Fisher zero-point has to be corrected for by either imposing oc H(z)~ 1 
or e g f oc H{z). We consider the former option unlikely, since this predicts spirals at z ~ 1.0 to be 
much bluer than observed (see Vogt et al. 1996). We therefore consider a scenario in which the 
galaxy formation efficiency is redshift dependent. Combining the observed Tully-Fisher relation 
with equation (^) yields 

w = 0.678 (-) -±1, (29) 

where L* = 1.0 x 10 10 h~ 2 L & . Using the Tully-Fisher relation of Giovanelli et al. (1997), converted 
to the F-band using V — / = 1.0 as an average for spirals (de Jong 1996b), yields 

e gf = 1.562 h 2 n (30) 

(i.e., Giovanelli et al. find 7 ~ 3.0). Here we assumed that the Hl-linewidth equals twice the 
circular velocity V200- 

The important factor here is 

H(z) 



ft A + (l_0 A -ft )(l + z)2 +i 2o(l + 2)3, (31) 

flu 
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which shows that the galaxy formation efficiency depends on the specific cosmology and the 
formation redshift. Here J7a = A/(3Hq), with A the cosmological constant. In the SCDM model 
with h = 0.5 and f2 = 1.0, equation ( ]30j ) yields e g f = 0.39(1 + z) 3 / 2 . Galaxies that were assembled 
at higher redshift thus formed more efficiently. More efficiently here means that a larger fraction 
of the baryonic matter available in the dark halo actually made it into the disk-bulge system. 
At z = this fraction is approximately forty percent, and this increases to hundred percent at 
z > 0.87. Halos of a given mass are denser at higher redshift, and higher density implies higher 
escape velocity and higher cooling efficiency, both of which could give a plausible explanation for 
an increase of the galaxy formation efficiency with increasing redshift. 

In what follows we consider two extreme cases of this galaxy formation (in) efficiency. In the 
first scenario, which we call the 'cooling'-scenario, we assume that only the inner fraction e g f of 
the baryonic mass is able to cool and form the disk-bulge system: the outer parts of the halo, 
where the density is lowest, but which contain the largest fraction of the total angular momentum, 
never gets to cool. In the second scenario, referred to hereafter as the 'feedback'-scenario, we 
assume that the feedback from star formation and evolution is such that each baryon in the dark 
halo, independent of its specific angular momentum, has a probability of e g f to ultimately end up 
as a constituent of the disk or bulge. The choice of either of these scenarios enters our equations 
in the function f(t,e g f) (see the Appendix). 



2.4. Stability 

As is well known, disks that are self-gravitating tend to be unstable against bar formation. 
Halos with smaller A yield disks with smaller scale- lengths (see equation ^0|). For a given 
disk-to-bulge ratio, there is thus a critical A such that for A < A cr it the disk will be bar unstable. 
Christodoulou, Shlosman & Tohline (1995) derived a stability criterion for bar formation based 
on the angular momentum content of the disk, rather than the energy content (such as the well 
known Ostriker & Peebles (1973) criterion). For a disk-bulge-halo system this stability criterion 
can be written such that disks are stable against bar formation if 

a = t^t < acrit, (32) 

with Vd and V c the characteristic circular velocities of the disk and the composite disk-bulge-halo 
system (equation ||J), respectively. This criterion is similar to the one proposed by Efstathiou, 
Lake & Negroponte (1982). As characteristic velocities we consider the circular velocities at 
R = 3Rd, since Mo et al. (1998) have shown that this is the radius where typically the rotation 
curve of the composite system reaches a maximum. For a thin exponential disk, approximately 80 
percent of the disk mass is located inside 3Rd, and 



Vi(3R d ) « 0.359 



(33) 
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Substituting equation (^) in (32), and using equation (20) we find that disks are stable if the spin 
parameter of their halos obeys 



with 

v^359 

^crit — • i-aOJ 

2a C rit 

This reduces to e cr it = 1.15 for stellar disks and e cr i t = 0.86 for gaseous disks, using a cr i t = 0.26 
and a cr i t = 0.35, respectively (see Christodoulou et al. 1995). This is in excellent agreement with 
the empirical value of e cr it = 1-1 found by Efstathiou et al. (1982) based on JV-body experiments. 
Throughout we take e cr it = 1.0 to take into account that disks are often made up of a mix of gas 
and stars. 



The critical value of A given by equation (34) depends on the galaxy formation efficiency 



e g f, and thus on the specific cosmology and formation redshift according to equation fl3"ip. 
Furthermore, since A cr it oc jT 1 , and jd is different for the cooling and feedback scenarios, A cr it 
depends on the particular model we choose for the galaxy formation (in)efficiency. In Figure 1 
this dependence is shown for three different cosmologies: SCDM, OCDM, and ACDM. The latter 
is a model with a non-zero cosmological constant with £Iq = 0.3, Q\ = 0.7, h = 0.6 and as = 0.97. 
Results are shown for both the feedback scenario (left panel) and the cooling scenario (right panel) . 
Figure 1 plots the logarithm of the critical spin parameter versus redshift. The dotted lines outline 
the probability distribution of the spin parameter (equation 0); the indices (indicated in the right 
panel only) denote the probability of finding a halo below that line. At low formation redshifts 
A cr it increases with redshift, due to the increasing galaxy formation efficiency with redshift. The 
turnovers to a flat line at higher redshift are caused by the requirement that e g f < 1.0. Note that 
in a SCDM cosmology with Zf OTm > 0.87 (for which e g f = 1.0) about 50 percent of all halos yield 
unstable disks, and thus require the formation of a bulge component in order to stabilize the disk. 
This increases to ~ 90 percent in models with 0,q = 0.3. Furthermore, it is apparent that A cr it is 
much less redshift dependent in the case of the cooling scenario, than in the case of the feedback 
scenario. The implications of this are discussed in Section 4. As can be seen, the results are very 
similar for the two low models. Since A cr it is the important parameter for what follows, we only 
focus on a comparison of SCDM with OCDM cosmologies, and note that results for the ACDM 
model are similar as for the OCDM model. 



2.5. Isothermal halo profiles 

It has recently been argued that the results of Navarro, Frenk & White concerning the 
universal density profiles of dark halos may be incorrect. Very high resolution simulations by 
Moore et al. (1997) indicate that the central density cusp of dark halos may be even steeper than 
r _1 . We therefore investigate the sensitivity of our results to the assumption of a NFW halo profile 
by comparing our results with those obtained considering isothermal halo profiles, i.e., halos with 
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a density distribution given by 

For such halo profiles the circular velocity, V c , is independent of radius and equal to the circular 
velocity at the virial radius, V2oo- The results for the isothermal density profiles can be obtained 
from the equations in this paper by taking f c = 1. We have rerun all models with an isothermal 
halo, and compared the results with those of the NFW halo profile. None of the results presented 
in this paper are significantly influenced by the assumption of a NFW halo profile. 



3. Towards understanding the Hq—R<i diagram 

One of the strongest suggestions for a smooth transition from spirals and SOs towards disky 
ellipticals along a sequence of decreasing disk-to-bulge ratio is provided by the Ho-Rd diagram, 
in which the central surface brightness of disks, no, is plotted against the disk's scale-length R<j 
(see Scorza &; Bender 1995). Scorza &; van den Bosch (1998) investigated the continuity of disk 
properties of a variety of disk galaxies and extended the diagram towards much higher fj,o and 
smaller Rd by including the nuclear disks discovered from Hubble Space Telescope photometry 
in a number of early-type galaxies (van den Bosch et al. 1994; Kormendy et al. 1996; van den 
Bosch, Jaffe & van der Marel 1997). This extended Ho~Rd diagram, already published in Scorza 
& van den Bosch (1998), is shown in Figure 2. It plots the central surface brightness (in F-band) 
as a function of the logarithm of the scale- length of a vast variety of disks, ranging from the 
extraordinary large disk of Malin I (data taken from Bothun et al. 1987), to the very small nuclear 
disks (data taken from Scorza & van den Bosch 1998). In addition, Figure 2 plots the parameters 
of embedded disks (solid circles, data taken from Scorza & Bender 1995 and Scorza et al. 1998), 
and of a combined sample of SOs (stars), HSB disks (open circles), and LSB disks (triangles). Data 
are taken from Kent (1985), de Jong (1996a), Sprayberry et al. (1995), de Blok, van der Hulst & 
Bothun (1995), and McGaugh Sz Bothun (1994). Figure 2 combines the parameters of disks that 
span 4 orders of magnitude in both scale-length and central surface brightness! 

Although there is considerable 'scatter', the /j,o(V)-Rd diagram of Figure 2 seems to reveal 
an almost linear relation. It is tempting to believe that this may teach us something about the 
formation of galactic disks. However, it should be noted that selection effects play an important 
role at the lower left part of the diagram: Most of the spirals and SOs in Figure 2 are taken from 
large diameter limited surveys, and are thus strongly biased towards the largest scale-lengths. 
Furthermore, one is always strongly biased against low surface brightness disks, especially if they 
are embedded in a much brighter spheroid. Even more, the fact that selection effects are different 
for the various data sets combined in Figure 2 prohibits a statistical conclusion on this lower left 
part of the diagram (e.g., see detailed discussions in de Jong 1996b and Dalcanton et al. 1997). 
However, the absence of disks in the upper right part of the diagram is real, and should in principle 
be related to the formation and evolution of disk systems. 

We now use our disk-bulge-halo models to investigate the locations of the different sorts of 
disks in the ^{y)-Rd diagram. First we note that spirals and SOs with the highest central surface 
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brightness for given scale-length lie approximately along a line of constant disk luminosity. Part of 
the upper boundary in the fj,o(V)-Rd diagram thus seems related to a cut-off in disk luminosity. 
This is to be expected from the cut-off in the galaxy luminosity function, which is well described by 
the Schechter (1976) function. This tells us that there are only very few galaxies with luminosities 
exceeding 3L* . We thus assume that the upper boundary is a line along which the sum of bulge 
and disk luminosities is equal to 3L*. The dotted line in Figure 2 (which is partly drawn as a 
solid thick line) corresponds to = 3L*. The position of a disk along a line of constant L<i is 
determined by the spin parameter A: disks originating from halos with smaller A have smaller 
scale-lengths, and higher central surface brightness. As shown in Section 2.4, there is a critical spin 
parameter such that for A < A cr it disks are unstable. The disk for which = 3L* and A = A cr it is 
indicated by a thick solid dot in Figure 2. If we want to build disks with scale- lengths smaller than 
this, we need a bulge to help stabilize the disk. The curved thick solid line, that originates at the 
thick solid dot and extends towards smaller scale-lengths, corresponds to a line with L^ + L^ = 3L*, 
A = A cr it, and with a disk-to-bulge ratio t running from t = oo (big solid dot) towards t = 0. The 
dashed line is the line with A = A crr t, t = oo, and running from 3L* (big solid dot) to = 0. 
Using these limits on the expected scale-lengths and central surface brightnesses of disks in our 
disk-bulge-halo models, we can distinguish 4 different regions in the no(V)-Rd diagram, labeled I 
to IV in Figure 2. Region I is the region where the disk luminosity exceeds 3L*. Region II is the 
region where < 3L* and A > A cr it; disks in this region are stable against bar formation even 
without a bulge. Disks in region III, however, require a bulge to be stable. Region IV, finally, is 
the region where systems with + L\ } < 3L* have an unstable disk despite the presence of the 
bulge. The position of the big solid dot, and therewith of the boundary between regions II and 
III, depends strongly on the formation redshift and cosmology (cf. Figure 1). Figure 2 is based on 
a SCDM cosmology with zero formation redshift. Although quite a number of assumptions have 
gone into our disk-bulge-halo models, the upper boundary (indicated by the thick solid line) nicely 
separates the region without disks, from the region most densely occupied by disks. Of course, the 
fact that the number density is highest just below the boundary is an observational bias effect as 
discussed above. One has to keep in mind that the (iQ—Rd diagram can be rather confusing since 
there are three parameters that determine the position of a disk in this diagram, namely the disk 
luminosity, the spin parameter, and the disk-to-bulge ratio. In the following we therefore focus on 
another diagram in which the luminosity dependence has been removed. 



4. Clues to the formation of bulge-disk systems 

We now ask ourselves whether we can learn something about the formation of disks and bulges 
from a comparison of our disk-bulge-halo models with real galaxies. We do this by calculating, 
for each galaxy in our combined sample, the spin parameter A of the dark halo which, under a 
given set of assumptions, yields the observed disk properties (scale-length and central surface 
brightness) . We thus use our formation scenario to link the disk parameters to those of the dark 
halo, and use the statistical properties of dark halos to discriminate between different cosmogonies. 
We start by choosing a cosmology and a formation redshift and we read in the disk-to-bulge ratio, 
t, and the scale-length and central surface brightness of the disk. From this we determine the 
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luminosities of the disk and bulge, which, upon using equations (|3|), (||), and (|30|), yield the total 
mass M and the particular density distribution of the disk-bulge-halo system. We then calculate 
A from equation (20) and plot this against the disk-to-bulge ratio t. This A is the value of the 



spin parameter that the halo would have had if the observed disk-bulge system formed out of that 
halo in the way envisioned here. One of the advantages of this diagram is that the dependence 
on luminosity is eliminated: the position of a disk in this diagram is virtually independent of the 
disk's total luminosity. 

In Figures 3 (SCDM) and 4 (OCDM) we plot the location of disks of LSB spirals (triangles), 
HSB spirals (open circles), SOs (stars), and disky ellipticals (solid circles) in a log [A] vs. log[i] 
diagram. The five small-dashed lines in Figures 3 and 4 correspond to the value of A below which 
we expect 1, 10, 50, 90, and 99 percent of the halos for a distribution function of spin parameters 
given by equation (0) with A = 0.05 and a\ = 0.7 (see also Figure 1). These lines are drawn to 
guide the eye. If our assumptions for the formation of disk-bulge-halo systems are correct, the 
distribution of inferred halo spin parameters should be similar to this log normal distribution. 

The thick solid lines in Figures 3 and 4 correspond to A = A cr i t (equation |$4|). Disks that lie 
below this line are unstable against bar formation. This critical value of the spin parameter has 
been calculated for a halo with a total mass of 10 13 M , but is virtually identical for other halo 
masses. The crowding of disks at t = 10 is due to the fact that we have assigned a disk-to-bulge 
ratio of ten to systems for which no value of t was determined (these are mainly a number of 
the LSB spirals). The panels on the left correspond to Zf orm = 0.0, the panels in the middle to 
Zform = 1-0, and the panels on the right to Zf orm = 3.0. The upper three panels correspond to 
models in which e g f = 1.0, i.e., in which all the baryons in the dark halo make it into the galaxy. 
The panels in the middle correspond to the cooling scenario, and the lower panels to the feedback 
scenario. The main difference between the upper and the middle and lower panels is that the 
former have A cr it independent of formation redshift (there is a non-zero but negligible dependence 
on Zf orm left due to the fact that the concentration parameter c of the NFW halo profiles depends 
on formation redshift), whereas A cr it increases with Zf orm in the middle and lower panels (according 
to the curves in Figure 1). 



Since Rj oc T2qq (equation pOfl) and r2oo = Q-^-V-zqq/ H{z), disks forming in a halo with mass 
M and spin parameter A are smaller for higher Zf orm . A higher formation redshift thus implies a 
larger spin parameter in order to yield the observed value of R^, something which is evident from 
Figures 3 and 4. Furthermore, disks seem to follow the A = A cr j t curve, i.e., disk-bulge systems 
do not have bulges that are significantly more massive than required by disk-stability. This is at 
least consistent with the inside-out bulge formation scenario proposed in Section 2.1, in which the 
formation of disk and bulge is self-regulating and governed by disk stability. 

The OCDM model with constant and near to unity galaxy formation efficiency (upper panels 
in Figure 4) can be ruled out. In this scenario systems with a large disk-to-bulge ratio, such as 
spirals, are rare; only ~ 10 percent of the dark halos is expected to yield such systems. The 
stability line A = A cr ;t crosses the average A of the spin parameter distribution at D/B ~ 1.0, 
and thus half of the halos is expected to form systems with a disk-to-bulge ratio less than unity. 
This is in conflict with the fact that the major fraction of field galaxies are spirals. Furthermore, 
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stable spirals can only form at high redshift (zf orm ;> 3), since for lower formation redshifts the 
spirals fall below the stability line. In the OCDM cooling scenario, spirals that formed recently 
tend to be slightly more stable, but still, the probability that a certain halo yields a system with 
a large disk-to-bulge ratio is rather small, rendering this scenario improbable. However, for the 
feedback OCDM cosmogony, a promising scenario unfolds: proto-galaxies that collapse at high 
redshifts are preferentially found in overdense regions such as clusters. The densities of these 
virialized halos are high, and this (somehow) causes the galaxy formation to be very efficient; 
nearly all the baryons in the dark halos make it into the galaxy. As a consequence, the systems 
that form require relatively small disk-to-bulge ratios (e.g., similar to SOs) in order to be stable. 
As the Universe evolves, and proto-galaxies with smaller over-densities (which are predominantly 
located in the field) start to turn around and collapse, the galaxy formation efficiency decreases, 
and systems with larger disk-to-bulge ratios are allowed to form. This scenario thus automatically 
yields a morphology-density relation, in which systems with smaller disk-to-bulge ratios (e.g., SOs) 
are preferentially to be found in overdense regions. Unfortunately, the disky ellipticals seem to not 
fit in this scenario. Even at formation redshifts of three, the average disky elliptical lies near the 
line with p(< A) = 0.01. 

The scenario outlined above for our OCDM cosmology also applies to the SCDM model with 
S7o = 1, albeit somewhat less restricting. After all, A cr i t only gets as large as ~ A. Consequently, in 
a SCDM Universe our model predicts less of a morphology-density relation than in a low density 
Universe; i.e., the magnitude of the morphology-density relation depends on the baryon fraction 

^bar/^0- 



4.1. Formation redshifts 

We define the typical formation redshift of a subset of galaxies as that redshift for which the 
average disk yields A = A. In Figure 5 we plot histograms of inferred formation redshifts for spirals 
(HSB and LSB disks combined), SOs, and disky ellipticals. The three dotted lines are drawn 
to guide the eye and correspond to redshifts of 1, 3, and 10. The inferred formation redshift is 
determined by calculating the redshift at which a halo with A = A = 0.05 yields the observed disk 
parameters. If a disk lies already above the A = 0.05 line at z = 0, a formation redshift of zero is 
assigned. The arrows in Figure 5 indicate the mean of the inferred formation redshifts for each 
sample. As can be seen, spiral disks are consistent with having formed fairly recently: (zf orm ) « 0.3 
with a mean formation redshift of zero. This conclusion was already reached by Mo et al. (1998) 
based on the same disk formation scenario as used here, but on somewhat different arguments. 
The distribution of Zf OTm is much more spread out for SOs, with a mean formation redshift of ~ 7. 
Note that misclassifications may be partially responsible for this large spread. Finally, the inferred 
formation redshifts for the disky ellipticals are unrealistically high, with a mean of ~ 38. We thus 
conclude that unless there is a significant loss in specific angular momentum of the material that 
forms these systems, disky ellipticals seem to not form a continuous sequence with SOs and spirals, 
at least not in the sense of the formation scenario investigated here. 

The histograms in Figure 5 correspond to the feedback SCDM cosmogony. However, the 
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results are very similar for the other cosmogonies discussed in this paper: although the details of 
the histograms are somewhat different, the averages are the same within the errorbars. 



As discussed above it seems unlikely that disky ellipticals can be fit in the simple bulge-disk 
formation scenario envisioned here. The small disk-to-bulge ratios in these systems renders an 
inside-out formation scenario for the bulges in disky ellipticals unlikely Two alternative formation 
scenarios for disky ellipticals have been proposed: (i) disky ellipticals are the outcome of gas-rich 
mergers (e.g., Bender, Burstein & Faber 1992; Scorza & Bender 1996), and (ii) embedded disks 
are formed out of internal mass loss from the old stars in the spheroid (Scorza 1993; Brighenti & 
Mathews 1997). 

Hernquist & Barnes (1991) have shown how in gas-rich mergers the gas becomes segregated 
in two components. The first component looses most of its specific angular momentum and 
collapses to form a dense central cloud. When converted into stars, these compact blobs of gas 
may well account for the high central densities in disky ellipticals (Barnes 1998). The second 
component gets propelled to larger radii in the tidal tails originating from the merger. This gas 
may either escape or reaccrete (through cooling) onto the newly formed 'elliptical' and form a new 
disk component, resembling the embedded disks. The numerous physical processes at play and 
the large parameter space available make it difficult to come up with detailed predictions for the 
resulting merger remnant. However, at present there seem to be no clear observations that rule 
out gas-rich mergers as the origin for disky ellipticals (but see Kissler-Patig, Forbes & Minniti 



In the internal mass loss scenario one expects a clear correlation between the mass of the 
disk and the mass of the spheroid. This correlation is indeed present and Scorza & Bender (1998) 
showed that typical mass loss rates for the old spheroid population are in excellent agreement 
with the typical ratio between the masses of the embedded disk and the spheroid. Another 
prediction from this scenario is that the spheroid and disk rotate in the same direction and that 
their specific angular momenta are similar (at least if we assume that the gas blown into the 
interstellar medium by the spheroid does not loose its specific angular momentum when cooling 
to form the embedded disk). Scorza & Bender (1995) found this to be the case from kinematic 
decompositions of line-of-sight velocity profiles of a number of disky ellipticals from which they 
determined estimates of the specific angular momenta of the embedded disk and the spheroid. 
Their data, which we plot in Figure 6, yields an average ratio of 



A third prediction associated with this scenario is that the disks are younger than the spheroids 
Recently, de Jong & Davies (1997) reported that disky ellipticals have higher H/3 line indices 
than the more luminous boxy ellipticals, and suggested this to be due to a young population, 
most likely the embedded disk (see also Brighenti & Mathews 1997). A final prediction for the 



5. Alternative formation scenarios for disky ellipticals 



1998). 




(37) 
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internal mass loss scenario is that the disks have surface brightness profiles that are steeper than 
exponential, since the gas does not cool from a uniform sphere in solid body rotation, but from a 
differentially rotating, centrally concentrated gas distribution (Scorza 1993; Brighenti & Mathews 
1997). Indeed, Scorza & Bender (1995) and Scorza et al. (1998) have found that the embedded 
disks often have surface brightness profiles that are significantly steeper than exponential. 

The masses, specific angular momenta, surface brightness profiles, and ages of embedded disks 
are thus all consistent with them having formed out of mass loss from the old stellar population of 
the spheroid. More accurate determinations of the specific angular momenta of embedded disks 
and their spheroids in a larger sample of disky ellipticals will proof helpful in constraining this 
further. 



6. Conclusions &: Discussion 

Understanding galaxy formation is intimately linked with understanding the origin of the 
Hubble sequence. An important clue is provided by comprehending the formation of disky 
ellipticals and SOs, simply because these systems form the transitional class from the classical 
ellipticals to the spirals. Disky ellipticals seem to form a continuous sequence with spirals and SOs 
towards smaller disk-to-bulge ratio, and it is thus tempting to believe that disky ellipticals, SOs 
and spirals all formed in a similar fashion. The aim of this paper has been to investigate a simple 
formation scenario for disk-bulge-halo systems, and to search for the main parameters and/or 
processes that determine the disk-to-bulge ratio and thus explain to a large extent the origin of 
(the major part of) the Hubble sequence. 

We considered the disk formation scenario originally proposed by Fall & Efstathiou (1980), 
in which the size of the disk, which is formed by cooling of the gas in a dark halo, is determined 
by the amount of angular momentum of the halo. We have extended upon this formation scenario 
by including a simple picture for the inside-out formation of an additional bulge component out 
of the inner, low angular momentum material. Stability arguments are used to suggest that the 
formation of the bulge is a self-regulating process in which the bulge grows until it is massive 
enough to allow the remaining gas to form a stable disk component. We do not describe the bulge 
formation in any detail but merely use empirical relations which allow us to describe the bulge by 
a single parameter, namely its mass. 

Each dark halo contains a fraction fibar/^o of baryons. We introduced a galaxy formation 
efficiency e g f which describes the fraction of those baryons that actually build up the disk-bulge 
system. The theory of spherical collapse coupled with the definition of a virialized halo predicts 
a Tully-Fisher relation of the form L oc as observed, with a zero-point that depends on the 
Hubble constant H(z). Recent observations, however, suggest that the Tully-Fisher zero-point is 
independent of redshift, implying that the galaxy formation efficiency is proportional to H{z). A 
physical explanation for this redshift dependence may be the higher escape velocities and cooling 
efficiencies at higher redshifts. 

For a combined sample of ~ 200 galaxies, varying from spirals to disky ellipticals, we 
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calculated the value of the halo's spin parameter which yields the observed disk properties under 
the assumption that disk-bulge systems form in the way envisioned here. We compared two 
cosmologies (SCDM vs. OCDM) and investigated the differences between assuming a galaxy 
formation efficiency of unity and two (extreme) scenarios in which e g f oc H(z): a cooling scenario, 
in which we assume that only the inner fraction e g f of the available baryons cools to form 
the disk and bulge, and a feedback scenario, in which each baryon, independent of its specific 
angular momentum, has a probability e g f of ultimately ending up in the disk or bulge. Our main 
conclusions are the following: 

• Disk-bulge systems do not have bulges that are significantly more massive than required 
by stability of the disk component. This suggests a coupling between the formation of disk 
and bulge, and is consistent with the self-regulating, inside-out bulge formation scenario 
proposed here. 

• If we live in a low-density Universe (Qq <^ 0.3), the only efficient way to make spiral galaxies 
is by assuring that only a relatively small fraction of the available baryons makes it into the 
galaxy, and furthermore that the angular momentum distribution of those baryons is similar 
to that of the entire system; i.e., the probability that a certain baryon becomes a constituent 
of the final galaxy has to be independent of its specific angular momentum. In the cooling 
scenario, most of the angular momentum of the system remains in the outer layers, and most 
halos form disk-systems with massive bulges, such as SOs. If, however, the galaxy formation 
efficiency is regulated as described by our 'feedback' model, a promising scenario unfolds: 
At formation redshifts ;> 3.0 the galaxy formation efficiency is unity, and systems that form 
build up a large bulge to support the disk that assembles around them. These galaxies 
resemble SOs. Galaxies that form later, at z ~ 0, no longer require a massive bulge, and 
spirals are preferentially formed. Coupled with the notion that density perturbations that 
collapse early are preferentially found in high density environments such as clusters, this 
scenario automatically predicts a morphology-density relation in which SOs are most likely 
to be found in clusters. In a SCDM Universe a similar, albeit less restrictive, mechanism is 
at work, which predicts a morphology-density relation of smaller amplitude. 

• A reasonable variation in formation redshift and halo angular momentum can yield 
approximately one order of magnitude variation in disk-to-bulge ratio, and our simple 
formation scenario can account for both spirals and SOs. However, disky ellipticals have 
too large bulges and too small disks to be incorporated in this scenario. Apparently, their 
formation and/or evolution has seen some processes that caused the baryons to loose a 
significant amount of their angular momentum. Merging and galaxy harassment (Moore et 
al. 1996) are likely to play a major role for these systems. 

Finally we wish to emphasize of few of the major shortcomings of the oversimplified formation 
scenario discussed here. First of all, we have neglected the fact that the merging of halos is an 
ongoing process, and that this is very effective in destroying disks (e.g., Toth & Ostriker 1992). 
Gas-dynamical simulations that do not involve the energy and momentum feedback to the gas 
from supernova explosions, stellar winds, UV radiation etc. produce galactic disks that are some 
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two orders of magnitude smaller than the observed spiral disks (e.g., Navarro & Steinmetz 1997): 
merging (and also harassment) are very effective in transporting angular momentum out into the 
halo, thus yielding more and more compact galaxies. This problem with forming galactic disks of 
proper dimensions is often referred to as the angular momentum problem. Despite our ignorance 
regarding the effects of merging, harassment, feedback, and (re)-ionization of the Universe, 
the observed sizes of (spiral) disks clearly suggest however, that the combine effect of all these 
processes is apparently such that the material that forms galactic disks has not lost much of its 
original angular momentum acquired from cosmological torques. The use of the Fall & Efstathiou 
disk formation scenario thus seems justified, despite the aforementioned shortcomings. 

In the past semi-analytical simulations of galaxy formation have been mainly based on the 
assumption that all bulges result from the merging of disk galaxies (e.g., Kauffmann, White & 
Guiderdoni 1993; Cole et al. 1994; Baugh, Cole & Frenk 1996; Somerville & Primack 1998). In 
this paper we have examined an inside-out bulge formation scenario, which should be regarded 
as complimentary to this merging scenario. Our formation scenario can account for both spirals 
and SOs, but fails to build systems that are even more bulge dominated. Although it has been 
shown that bar-instabilities can lead to the formation of small (and close to exponential) bulges, 
a more detailed study is required to investigate whether the inside-out bulge formation scenario 
discussed here can yield more massive bulges as well. It is at least intriguing that the critical spin 
parameter is of the same order of magnitude as the typical spin parameter for halos, suggesting 
that a significant fraction of halos will yield unstable disks, unless part of the baryonic material 
is transformed into a bulge component as suggested here. So despite the clearly oversimplified 
nature of the formation scenario envisioned here, it may provide a useful framework for future 
investigations of galaxy formation. 
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A. Angular momentum distribution of the disk 



Mestel (1963) noticed that disks of spirals have specific angular momentum distributions 
similar to that of a uniform sphere in solid body rotation. Exponential disks are thus a natural 
consequence of the Fall & Efstathiou disk formation scenario, if we assume that proto-galaxies 
have a specific angular momentum distribution similar to that of a uniformly rotating, uniform 
sphere, and that the material that forms the disk does not loose its specific angular momentum 
(e.g., Gunn 1982; Dalcanton et al. 1997). 



The total angular momentum of a uniformly rotating, uniform sphere with radius a and 
constant angular velocity u> is 



J = -Mcoa, 
5 

and the mass with specific angular momentum less than j is given by 



M(< j) = M 
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Since the baryonic and dark matter are well mixed, the angular momentum distribution for 
the baryonic matter component is identical to that of the total mass of the system, and 
^bar = (^bar/^oM- The bulge forms out of the low angular momentum material with a cylindrical 
radius R < -Rbuigc- Solving Mb ar (< j) = Mj, yields 
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and the total angular momentum of the bulge is 
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which is the total angular momentum of the intersection of a solid sphere with radius a and a 
solid cylinder with radius -Rtmige < a (where the sphere and cylinder have a common center). 
The disk forms out of the material with -Rtmige < R < a, and has a total angular momentum 
Jd = <^bar — Jb = jd,J ■ We assume that the specific angular momentum of this material is conserved 
during disk formation. 

In the feedback scenario, in which each baryon has a probability e g f of ending up in the final 
disk-bulge system, we find upon substituting -Rbuige (equation [A3]) in equation (|A4|) , and with 
the use of equations (||) and (§]) 



Jd 



[m d + m b ] 



l + t 



l + t 



2/3' 



f(t,e gf ) [m d + M b ] 



(A5) 



In the cooling scenario, only the baryonic material with r < r coo i cools to form the disk and 
bulge, where 



^cooi = e„( 3 a, (A6) 
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and one finds 



3 d = e 2 J?{m d + m h )- t 



+ t 



t \ 2 / 3 
1 + J 



= f(t,e gf )[m d + M b ]. 



(A7) 
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Fig. 1. — The critical value A cr ;t of the halo's spin parameter (plotted logarithmically) as a function 
of redshift z for both the feedback scenario (left panel) and the cooling scenario (right panel). 
Results are plotted for three different cosmologies: SCDM (solid line), OCDM (dot-dashed line), 
and ACDM (dashed line). The five dotted lines are labeled by the probability of finding a halo 
below that line for the spin parameter distribution discussed in Section 2.2. The increase of A cr it 
with increasing redshift is due to the increasing galaxy formation efficiency. When this efficiency 
becomes unity, A cr it becomes constant with redshift as apparent in the figure. Note the similarity 
between the OCDM and ACDM cosmologies compared to the SCDM scenario. 

Fig. 2. — The Ho-Rd diagram for a large variety of disks: LSB spirals (triangles), HSB spirals 
(open circles), SOs (stars), disky ellipticals (solid circles), and nuclear disks (squares). The dotted 
line (partly plotted as a thick solid line) corresponds to the line with = 3L* . The dashed line 
corresponds to A = A cr i t for the case without bulge (i.e., t = oo). The curved thick solid line 
corresponds to A = A cr it and + = 3L*, with t running from zero to infinity. These lines 
border 4 regions, the meaning of which are discussed in the text. 

Fig. 3. — Results for a SCDM cosmology. Plotted are the logarithm of the spin parameter versus 
the logarithm of the disk-to-bulge ratio {D/B = t). Solid circles correspond to disky ellipticals, 
stars to SOs, open circles to HSB spirals, and triangles to LSB spirals. The thick solid line is the 
stability margin A cr it ; halos below this line result in unstable disks. As can be seen, real disks avoid 
this region, but stay relatively close to the stability margin, in agreement with the self-regulating 
bulge formation scenario discussed in Section 2.1. The dashed curves correspond to the 1, 10, 50, 
90, and 99 percent levels of the cumulative distribution of the spin parameter according to equation 
(|7|) with A = 0.05 and a\ = 0.7. Upper panels have a galaxy formation efficiency equal to unity, 
middle panels correspond to the cooling scenario, and lower panels to the feedback scenario. Panels 
on the left correspond to Zf orm = 0, middle panels to Zf 0Tm = 1, and panels on the right to Zf orm = 3. 

Fig. 4. — Same as Figure 3, except now for the OCDM universe with Qq = 0.3. The main difference 
is that A cr ;t reaches higher values here. Both in the case with e g f = 1.0 (upper panels), and in 
the cooling scenario (middle panels), ~ 90 percent of the halos is expected to result in unstable 
disks unless a significant bulge is formed (in the cooling scenario this reduces to ~ 75 percent for 
-Zform = 0). The majority of the halos yields systems with a significant bulge such as SOs, and both 
these models can thus be ruled out, given the relatively large numbers of spirals in the present 
Universe. In the feedback scenario (lower panels) stable spiral systems can form, albeit only at low 
redshifts. At higher redshifts systems with smaller disk-to-bulge ratios form, such as SOs. 

Fig. 5. — Histograms of the (normalized) distribution of formation redshifts (plotted as log[l + 
Zform]) determined as described in the text for spirals (left panel, HSB and LSB spirals combined), 
SOs (middle panel), and disky ellipticals (right panel). The arrows indicate the means for each 
sample (the mean for the spirals is zero). The number n of galaxies in each of the separate samples 
is indicated in the corresponding boxes. The dotted lines indicate formation redshifts of z = 1, 
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z = 3, and z = 10. These histograms correspond to the SCDM model with feedback, but are similar 
for other cosmogonies. 

Fig. 6. — The specific angular momentum J/M (in kpc km s -1 ) of the disk versus that of the bulge. 
As can be seen, within the errorbars the specific angular momenta of disk and bulge are roughly- 
similar. (Data taken from Scorza & Bender 1995). 
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